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Abstract 



The Hamilton-Jacobi equation in the sense of Poincare, i.e. formulated in the 
extended phase space and including regularization, is revisited building canonical 
' transformations with the purpose of Hamiltonian reduction. We illustrate our ap- 

^ , proach dealing with orbital and attitude dynamics. Based on the use of Whittaker 

and Andoyer symplectic charts, for which all but one coordinates are cyclic in the 
Hamilton-Jacobi equation, we provide whole families of canonical transformations, 
among which one recognizes the familiar ones used in orbital and attitude dynamics. 
' In addition, new canonical transformations are demonstrated. 

(N 
^ . 

cn 1 Introduction 

^ : Transformation of variables, including or not the change of time used as evolution pa- 
Q\ [ rameter, is a basic tool in almost any problem of physics. A recent illustration of the 
I usefulness of this technique, which original application goes back to Laplace and Euler, 
■ may be found in [19]. Changes of variables belong to Geometry and, therefore, do not 
^ ■ depend on the physical system to which they may be applied. However, transformations 
are commonly derived for specific purposes and, hence, they are tightly related to the 
particular problem they are applied to. 

A renewed interest on transformation of variables tailored to specific problems has 
recently emerged in three areas. On one hand there is a continuous search for more efficient 
numerical integrators demanded in many applied fields, ranging from research groups 
endowed with powerful computing facilities, to research involved in very long periods of 
time (see [5l [16] and references therein). On the other hand the Hamilton-Jacobi (H-J) 
equation is receiving new attention in extended phase space, connecting this technique 
with the variational principles of Lagrangian and Hamiltonian Mechanics (see [31] and 
references therein). Finally a third area, which has its roots in Poincare's [30] fundamental 
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problem of dynamics, deals with Hamiltonian systems of the form Ti. = Ti.Q + Tii where 
Ho defines an integrable system, and Hi is considered a perturbation. The H-J equation 
is built from Ho by searching for a generating function such that the new canonical 
variables satisfy some requirements of perturbation theories (see, for instance, [37t \T7\). 
The content of this paper focuses in this last area. Finally note that the H-J approach 
has been used by many authors, along the history of this equation, looking for variables 
leading to separability. We do not consider that problem as such, although it might be 
included. 

There is a long tradition in Classical and Celestial Mechanics — nowadays a common 
ground in many other fields — of using both canonical transformations and changes of 
independent varible dubbed as regularization. The H-J equation is the usual approach, in 
which the canonical transformation is obtained from either a principal or a characteristic 
function |18j. 

The characteristic function approach does not rely on the vanishing of the new Hamil- 
tonian. It only requires a non-zero transformed Hamiltonian that is cyclic in all variables 
but depends on one or several momenta. The adequate selection of the new Hamiltonian 
will make the transformation successful for specific purposes. Thus, for instance, if one 
wants the transformed Hamiltonian to retain the topology of the original problem, prob- 
ably several momenta must be retained in the new Hamiltonian. On the contrary, useful 
transformations for perturbation theory may look for transformed Hamiltonians that de- 
pend only on one momenta, thus simplifying the subsequent application of perturbation 
methods (which, in turn, have limited application to specific solutions as, for instance, 
quasi-periodic motion). 

Sentences like "after several trials" , which are frequent in the literature, might produce 
the impression that succeeding in the choice of the reduced Hamiltonian is a matter of 
educated guess or alquimia. However, the new Hamiltonian does not need to be specified 
since the beginning and, quite on the contrary, it may be handled formally throughout the 
process of computing the transformation [35]. Therefore, Hamiltonian problems solved by 
the H-J equation give rise to whole families of canonical transformations from which the 
selection of a particular one requires only the materialization of the new HamiltoniaE0. 

We propose to obtain canonical transformations from a general form of a H-J equation 
that is formulated in the extended phase space. In principle, the only requirement imposed 
to the Hamiltonian function ([6]) is that Hq has to be integrable. More precisely, the 
two families studied are presented using symplectic charts where the unperturbed part 
defines a 1-DOF. In other words, we only deal with canonical transformations derived 
from integrable Hamiltonians that are cyclic in all but one variable, and perform two 
specific applications. In attitude dynamics, we find a family of transformations that 

^To make more clear our point from [35] (p. 412), where E denotes the Hamiltonian in the new 
variables, we quote "We are still free to choose the functional form of E. A convenient (and conventional) 
choice is E{pQ,pi,p2) = —m^'^/(2p'^). With this choice the momentum po has dimensions of angular 
momentum, and the conjugate coordinate is an angle." As we will show this is not the only way of 
reasoning. Indeed, we arrive to the same choice when we impose the function called Kepler equation as 
one of the expressions defining the transformation. In other words, when the mean anomaly is chosen as 
one of the new coordinates. As Delaunay showed, in doing so he avoided the presence of mixed secular 
terms in this theory of the Moon. For a recent study on these variables from a geometric point of view 
the reader should consult 
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affords for complete reduction of tlie Euler-Poinsot problem. In tliis case, we clieck tliat 
Sadov's [321 [33], Kinosliita's plj, Deprit's [I2], and tlie recent Fukusliima's reduced 
Hamiltonians are just particular cases of our family. Notably, our approach provides 
explicit transformations from Sadov's and Kinosliita's elements to Andoyer variables |3j|, 
thus overcoming the major objection to both sets of elements of being constructed through 
implicit relations. Besides, we show that other selections of the new Hamiltonian provide 
new canonical transformations that may help in simplifying perturbation algorithms in 
attitude dynamics. 

In orbital dynamics we obtain three different families of isochronal canonical trans- 
formations — derived from Henon's isochronal central potential [211 El — each one linked 
with a different form of the regularizing function. From our families we recover a variety 
of canonical transformations in the literature, ranging from historic ones like Delaunay 
[H], Levi-Civita [23 EE], or Hill [22], to the recent transformation of Yanguas [37]. In 
addition, we show that the set of Delaunay elements yields the simplest transformation 
of the isochronal family that reduces the Hamiltonian of the two body problem to one 
momentum and, besides, preserves the original form of the Kepler equation. 

The two charts used enjoy similar characteristics. Both express the unperturbed part 
of the Hamiltonian families as genuine 1-DOF systems due to the central and axial sym- 
metries that both possess. They belong to the category of superintegrable systems, but 
not maximally superintegrable (see Fasso [14] and Ortega and Ratiu ^29j); only for par- 
ticular choices of the parameters both systems fall into this last category. Moreover, from 
the computational point of view, these charts introduce great simplification in intermedi- 
ary expressions which, in turn, lead to a deep insight about each family. Many authors 
still rely on Eulerian variables when dealing with rotational dynamics (see for instance 
Marsden and Ratiu [28j) , and spherical variables for orbital type motions (see Goldstein 
et ai. [18], Jose and Saletan [23] and Sussman and Wisdom [35]). 

In spite of the fact that the two families considered are independent and consequently, 
they are studied in a totally selfcontained manner, we make both part of the paper because 
this allows to obtain a better insight in the way we approach the H-J formalism. A final 
comment is due. Although both families have parameters, we have made the analysis 
for the generic case only. With respect to the orbital case, as an illustration, we have 
added details when we restrict to Keplerian systems, showing the connection of the new 
variables with the anomalies. We would like to make clear that we are not doing a 
review of regularizations. In other words, it is left for the reader to try other canonical 
transformations within the general scheme presented here. 

2 A General Form of the Hamilton-Jacobi Equation 

We only deal with Hamiltonians of the type 

;C(a;o,a;,Xo,X;/i) = (Xo + 7^)x, (1) 

where coordinates and X = (Xi,...,X„) conjugate momenta; Xq 

is the independent variable and Xq its conjugate momenta in the extended phase space 
formulation, in which we restrict to the manifold /C = 0; /i is a vector of parameters, and 
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the Hamiltonian 7i as well as the "regularizing factor" x i^^iy depend on all or some of 
the parameters defining /z: 

X = x(a;o,a;,Xo,X;/x), H = n{xo, x, X; n). (2) 

A dash in the place of a variable is used to remark that the corresponding variable is not 
present. 

Hamilton equations are 



dxo 


dJC 


dXo 


die 


dr 


dXo' 


dr 


dxQ 


dxi _ 


die 


dX, 


ok: 


dr ~ 


dX,' 


dr 


dxi' 



(3) 

(2 = l,...,n) 

where r is the evolution parameter of the flow. Note that the first of the previous equations 
reads 

dxo , . 

d7 = ^ (4) 
which tells the function x ought to verify that x > in its domain. Moreover, in the case 
of conservative systems Ti = H{—,x, —,X), Xq is an integral, and the manifold = 
may be also seen as Xq = —Ti = constant. 

We are interested in canonical transformations 

(xo, X, Xo, X) (?/o, y, Yq, Y) (5) 

in the sense of Poincare. More precisely, we look for transformations such that they 
simplify Hamiltonian systems defined by functions Ti. which can be written as 

n = no + ni (6) 

where Ho = Ho(— , x, — , X) defines an integrable system, and Hi = Hi{xo, x, X) is a per- 
turbation. Specifically, we focus on canonical transformations such that the new Hamil- 
tonian K = Kq + Ki satisfies 

Ko = (Xo + Ho)x = <^>(-, Y), (7) 

i.e. the full reduction of the unperturbed part is carried out. This Eq. (JTj) is what we refers 
as the variant of Poincare to the H-J equation; the classical case chooses $(— , — , Iq) Y) = 
0. 

The transformations are defined by 

X. = ^, ,. = ^, . = 0,...,n (8) 

derived from a generator W = W(xo, x, Yq, Y, /i) that is a complete solution of the gener- 
alized H-J equation 



dW ^, ( dW' 
VHq\x, — — 

OXq \ OX 



x(^x,a:o,^,^j=$(lo,n (9) 
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Thus, the Hamihonian /C in the new variables will take the form 



K; = /Co + /Ci = $ + T^ix (10) 

where Tii and x ^i-re expressed in the new variables. Note that in what follows we take 
Hi = 0. 

In this paper we limit to generators of the form 

W= Yl x,Y, + n{xn,Yo,Y), (11) 

0<j<n-l 

and regularizing function x = xi^n, Xq, Xi, . . . , Xn-i)- Hence, from Eq. ([9]) we may write 
Hofx,yi,...,y. J^]- ^(^o>^i>---^^-) y^. (12) 

Depending on the form of Hq, Eq. (fT2l) may be solved for dTZ/dxn and, therefore, 71 
is computed from a quadrature, which solution will depend on the choices made for $ and 
X- Other possibilities are under investigation [15]. Note that, in fact, there is no reason 
why we should impose on Ho to be cyclic in xq. What we have presented above, properly 
adapted, remains valid if we lift that constraint. This is referred in the literature as 
nonautonomuos systems; the driven oscillator, the relativistic particle, etc are just simple 
examples within that category (for more recent systems of interest see Struckmeier [3T]). 
In the families we will study below the possible presence of Xq occurs in the perturbing 
part. 

Apart from the general case of Eq. f[T^ . transformations non-based on the homoge- 
neous formalism (% = 1, $ = Iq + ^), adopt the simpler formulation 

x„, Fi, . . . , ^ j = ^(Fi, ...,Yn). (13) 

In the two types of Hamiltonian systems we study below, instead of (xq, x, Xq, X) we 
will start each Section considering briefly the two symplectic charts on which our study 
relies. In doing so we will use the standard notations, although there are more than one: 
Andoyer variables [3] for rotational motion and Whittaker (polar-nodal) variables [36] for 
orbital dynamics. 

With respect to the new variables (t/q, 2/, ^o? Y), we will use (ci, 7, f , m, D, F, T, U) for 
the rotational families and (/, g, h, u, F, G, H, U) for the orbital ones. For the different 
sets within each family, we do not find necessary to distinguish among them, like using 
subindices, etc. Each one distinguishes itself by the expression taken by JCq. 



3 Rigid body transformations 

Within the frame set up by Eqs (I7l) and (IHl) we will consider first canonical transformations 
derived from free rigid body Hamiltonian function. As it is well known, meanwhile Euler 
of variables continues to be the ones used in classical and recent books (see [H], [28]) at the 
same time in research papers is already customary to consider those systems formulated 
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in Andoyer variables. More precisely, the Hamiltonian of the free rigid body takes the 
form 

n = -, u, -, M, N; tti) = i (ai sin^ u + a2 cos^ u) (M^ - iV^) + | ag N\ (14) 

where < 03 < 02 < are parameters. Then, according to our notation Eq. ([T]), we 
deal with the Hamiltonian 



T + i (ai sin^ u + a2 cos' z/) (M^ - A^') + 1 N^] x{v, A, M, T), (15) 



where (A, /i, z/. A, M, A^) are Andoyer variables [3], and M' > A^' making the Hamiltonian, 
(fT^ . strictly positive and Ix^^s^^ ^ 'H < Ix^H^^- As we said in what follows we 
consider only the generic case 03 < 02 < Oi. We do not need to give more details on those 
variables here. The interested reader may consult the classical note of Deprit ^ or in 
more geometric terms Fasso 



3.1 The general structure of the transformation 

We look for canonical transformations 

(A, /i, z/, t. A, M, AT, T) ^ (d, 7, V, n, D, T, T, ?7) (16) 

that convert ( fT5l) in a certain function $(D, F, T,f/) depending only on momenta. The 
transformation will be defined through a function S = S{X, fi,i>,t, D,r,T ,U) in mixed 
variables such that 

d = SD, 7 = '5r, v = Sr, u = Su, A = Sx, M = S^, N = S,, T = St. (17) 
where we use the notation Sx = dS/dx. Thus, from (fT5|) we set the H-J equation 

j (ai sin' u + a2 cos' z/) (5j - Sl) + | 03 Sl + 5^] x(^^, 5a, <S^, 5*), = $(/^, F, T, U) (18) 
and iS is chosen in separate variables as 

5 = f/t + TA + F/i + W(z/,D,F,T,[/) (19) 
because t. A, and u are cyclic in (ITS!) . Then, 

W = T r v/Q(z.,D,F,T,t/)dz/, (20) 

where 

ai sm z/ + 02 cos^ ^^ — 03 f y^ / 

and a/Q must be real for all z/; therefore, $ < x (^02 F^ + f/). 
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If we assume that the regularizing factor depends only on the coordinate z/, the trans- 
formation ( ITTl) are 

d = ^$dX3, (22) 

2U /2$ 1 \ 

7 = /i + Ji + Y^J2- (^^--$rjX3, (23) 

V = A + ^$tX3, (24) 

u = t-lj2 + i$f,X3 (25) 

N = tJq, (26) 



where 



T = U, A = T, M = r, (27) 
Ji = /'^ ./q dz., = r% dz., Xs = ^ dz/. (28) 



3.2 The case x = 1 

We only deal with transformations with x = 1- Therefore, X2 = X3. 

Note that a restriction of Andoyer variables is |A^| < M, which, because of ( !26|) . 
further constrains the values of the new Hamiltonian to be into the positive interval 
$ G (|a3 r^, r^). In addition, one should be aware that the sign of a/Q must be taken 
in accordance with the sign of N . 

Again, the transformation depends on the integration of two quadratures, which closed 
form solution requires well known changes of variables. Thus, introducing the parameter^ 

e=^i^, / = ^^^, (29) 

(1\ — O3 0,2 — C^S 

and the functions 5 = S{D, T,T,U; Oi, 02, a^), m = m{6), defined by 

ai — 02 



P ai -2($ -f/) 



P, (30) 



1-e T^ai-2{<^-U) 02-03 ' ^ ' 



the quadratures in (128!) are solved to give 



mF{ip\m) - {m + f)U{ip,-f\m) 

' ~ 7U ' ^ ^ 

X2 ^ "^^(^'-^ (33) 

v/a2-a3Vr2ai-2($-f/) 



^ For the sake of linking with recent hterature we used Fukushima's notation [T7]. Deprit [T^ calls 
6 = a^, £ = 0:3, and m—k^, and Sadov |32| names m = A and f — k^. 

7 



where F{il)\m) is the eUiptic integral of the first kind of modulus m which amplitude ip is 
unambiguously defined through 



cos u 



sin ip 



a/1 — e cos^ tjj 



sm u 



Vl — e cosip 
a/1 — ecos^ ip^ 



dz/ 



1 — e cos^ 



d^/^. (34) 



On the other hand, Il{ip, —f', fn) is the elliptic integral of the third kind of modulus m, 
characteristic — /, and amplitude ip. Alternatively, one can write 



;i-5)n(</),5|m). 



where the new amplitude (p is unambiguously defined through 
cos (p . ^/l — 6 sin (p 



cos u 



sm u 



1 — 5 sin^ (p 
If we further introduce 



dz/ 



6 sin^ 



r 



T^ai - 2($ - U) Va2 - ag «i - ^2 ^2 - ^3 



1 — 6 sin^ 



7' 



(35) 



d(p. (36) 



then, the transformation equations above, (I22l) - fl26|l . adopt the compact form 



d 

7 



^DF{ip,m) 
(r ag - $r) F{^, m) + r (ai - as) n(^, -/|m) 



r 

cr 

a r/ $ - f/ 

"-fK^— - 

A + ^<l>TF(^/',m) 



$rj F{^p,m) - ( 2 



$ - f/ 



ra2 n(0,5|m) 



1)F( 



t 



■ (i. 



D 



e{l — 6 cos^ z/) 
^ S (1 — e cos^ z/) 



'1 — m sin^ t/'). 
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(37) 



(38) 
(39) 
(40) 

(41) 

(42) 

(43) 



where we maintain two expressions for 7 to ease comparisons with transformations in 
the literature. Remark that $ remains arbitrary. Therefore, fl37p - (H5]) define a family of 
canonical transformations that provide complete reduction of the Euler-Poinsot problem. 
Note, besides, from P2l) that reduced Hamiltonians of the form $ = f/ + ^(-D, F, T) are 
required for transformations preserving the time scale. 

Despite the family of transformations above has been obtained from a specific dynam- 
ical system, the rigid-body problem, once it has been derived it may be applied to any 
problem we wish. However, successful transformations for different problems are closely 
dependent on the selection of $. Thus, for instance, if we are to apply one of the trans- 
formations of the family (!38|) - (l43|) above to reduce the rigid body dynamics or to study 
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a perturbed rigid body, we should consider specific facts. Because neitfier A nor T are 
altered by the transformation it seems natural to choose $ = t/ + \E'(-D, F). Then, u is 
the time, v, T, U, T and D, remain constant, and 

7 = 7o + $rM, d = do + ^DU. (44) 

Because the super-integrability of the Euler-Poinsot problem the closure of a generic 
trajectory is a two-torus |4|, Therefore, \l/ should depend both on D and F if we want 
to retain the topology of the problem. Specific selections of depending only on D 
— whose appearance is mandatory for equation (138!) to be defined — will produce just 
periodic solutions. Even particular selections \E' = \E'(D,F) may reduce the solutions to 
just periodic orbits. For instance, if we impose 

i^r=J ^D, (45) 

i, j, integer, any computed $ = t/ + "^{iD + j F) with arbitrary \E' will provide only 
periodic solutions. 

We discuss here some other possible choices of Without any application in mind, 
one criteria might be 'to simplify' the coefficients of the variables in the transformation. 
Thus, for instance, a = D in fl37|) implies 



FV \ 

$ = ai (46) 

2 V CL2- az) 

where F has dimensions of momentum and D of inverse of momentum of inertia. This 
Hamiltonian simplifies the transformation equations, (l38|) - fH3l) . to 

F ^2 

d = F(^,m) (47) 

a2 - as 

7 = fL + Dlai-a2 ]u(4>,S\m) (48) 

V a2 - as / 

V = \ (49) 
u = t, (50) 



N = TD 



m sin^ ip 



\ (ai - as) (a2 - as)' 
To discuss other simplifications, we find it convenient to introduce 



(51) 



^1 = ^' ^2= (2^— ^-^r):^, ^3 = 1^, (52) 



V F / $£) $D 

and consider the possible choices of $ that make Si {i = 1,2,3) independent of the 
momenta. 
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3.2.1 Case si = si o and S3 = S3 constant. 

Then, from the first and third of ( l52l) and from ( 1371) . we get 

■53,0 



and, therefore 



■Sl,0 



(53) 



(54) 



with \1/ arbitrary. 

Because the topology of the problem, we can add the condition $t = 0, and therefore 
53,0 = 0. In addition, if we want to preserve the time scale (H2l) . we restrict to the case 



$-f/ = ^r^-^^(r,/^)^ 



(55) 



Sadov's transformation. The choice 



^ = Si,o \ D 

V «2 - 03 



gives 



$ - t/ 



ai — 03 f D 



r 



ai 



ai 



1 - 



ai — 03 K 



that is Sadov's [32] reduced Hamiltonian in which = f and 



(56) 



1 . 



(57) 



Kinoshita's case. If we instead choose 



ai (ai + a2 — 203) 
\ (^2 - 03) (ai + 



we find 

ai + 02 



7^ 



1 / ai + 02 



{51 



as 



(ai + 02) < 0, (59) 



where Ti. is Kinoshita's reduced Hamiltonian [24, see Eq. (10') on p. 433]. Alternatively, 
the scaling may be avoided in our general formulation by choosing a regularizing factor 
X = 2ai/(ai + 02) independent of the variables. 

Thus, the major objection to Sadov's and Kinoshita's elements of being related to 
Andoyer canonical variables through implicit transformations is easily circumvented using 
the general transformation ( l371) -( l43|) above particularized for either Sadov's or Kinoshita's 
Hamiltonians, Eq. (!56|) and (!59|) respectively. 
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Fukushima proposal. If we choose now 



VI/ = si,o Vn - D\ (60) 

we get Fukushima' s proposal [17] 

^ = U + ^aiD\ (61) 

that, depending on less than the required momenta, constrain the topology of the trans- 
formed Hamiltonian to periodic solutions only. Nevertheless, this way of proceeding may 
be perfectly adequate for a perturbation theory, as shown in [37] . 



TR-type mapping. In a similar way to Fukushima, we propose to choose 



* = ^i,oJ-^(r-D), (62) 

leading to 

^ -U = aiD (r - Id) (63) 



that is formally equal to the reduced Hamiltonian of the TR-mapping [Ml [TT] given in 
Eq. (11251) below. It is quadratic, depends on two momenta — as it should be for the more 
general case of a super-integrable three degrees of freedom problem with four independent 
integrals [3] — and the momenta keep the dimensions of the original problem because the 
Hamiltonian is multiplied by the inverse of the momentum of inertia ai. 

Then, the reduced system has four elements {v, D, F, T) and two variables that evolve 
linearly with time 

7 = 7o + oi -D i, d = do + ai (T — D) t, 
the condition T/D rational providing the subset of periodic solutions. 

3.2.2 Case S2 = ^2,0 and S3 = S3 constant 
From the second and third of flS21) we get 

S2,o$D = — $T = 2^^-<l'r ^ $ = f/ + F2^(D-S2or + S3,oT,f/) (64) 

s3,o r 

where \Ef is arbitrary. 

A simple choice satisfying the condition in (1641) could be $ = f/ + F^ (Z) — S2,o T), for 
instance. Particular cases in the literature make S2,o = 0. Thus, Deprit and Elipe [12] 
choose 

1 F2 

<I>-U= , (65) 

Other solutions could be, for instance, 

^-U = T^D^, ^-U = T^D. (66) 
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4 Orbital Canonical transformations 



What we have said at the beginning of the previous Section for canonical variables in 
rotational dynamics, it appears again in orbital variables. Thus, meanwhile Hamiltonian 
particle dynamics still happens to be introduced in spherical canonical variables (see for 
instance [6], [18], [23], [35]), at the same time in space research another canonical set 
of variables is widely in use; we refer to nodal-polar variables, already considered for 
planetary theories almost a century ago [22j. The main difference of both set of variables 
is that nodal-polar variables carries out a double reduction (axial and central symmetries) 
of Kepler type systems, meanwhile spherical variables only deals with the axial symmetry. 

As in the search of orbital canonical transformations we build up our H-J equation 
based on nodal-polar variables, we briefly give a description about them. Taking a ref- 
erence frame (0,61,62,63), they are introduced in a natural way associated with the 
"instantaneous plane of motion" whose characteristic vector is x x X. The nodal-polar 
variables (r, i/, ^) are defined as follows: the variable r is ||x||; the angle v called the as- 
cending node, is defined by cos z/ = 61 ■ £ and sin u = e2 ■ i, where £ is the unit vector 
defined by 63 x (x x X); the variable ^ is a polar angle in the plane of motion giving 
the position of x reckoned from the vector i. The variables {R, N, B) are the correspond- 
ing conjugate momenta. For more details and explicit expressions of the transformation 
{x,y, z, X,Y, Z) — * (r, z/, 6', i?, iV, 9), called by some authors Whittaker transformation, 
see for instance [TO] . 



4.1 The general structure of the transformations 

Given the function 



n 



+ V(r)+T 



x(r,e,iV,T) 



(67) 



where (r, 6, u, t, R, B, A^, T) are Hill or polar-nodal variables in the extended phase space 
and V only depends on distance, we look for a canonical transformation 

(r, 9, u, t, R, e, N, T) ^ (/, g, h, u, F, G, H, U) 

that converts fl^TI) in a certain function $(-F, G, if, U) depending only on the momenta. 

The transformation will be defined through a generating function in mixed variables 
S = S{r, e, u, t, F, G, H, U) such that 

f = Sf, g = <Sg, h = Sh, u = Su, R = Sr, 6 = S0, N = S^, T = Sf (68) 

Then, from fl^ we set the H-J equation 



x{r,Se,S,,St) = ^F,G,H,U) 



(69) 



Because t, 6, and u are not present in ( 167|) . the generating function may be chosen in 
separate variables 

S = Ut + Hu + G9 + Wir,F,G,H,U) (70) 
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Then, 



1 (w^ + 91] + v(r) + U = ^(MJLin (71) 

2 V'^'^ r^J^'^^ '^"^ x{r,G,H,U) ^' ' 



that can be solved for W by a simple quadrature: 



where 



r ^Q{r,F,G,H,U)dr, (72) 



X (r, G, H, U) 
Therefore, the transformation is 

/ = $fX3, (74) 
g = 9 + GIi + <t>Gl3-<l> r ^^dr, (75) 

h = i. + $^X3-$ r^^dr, (76) 



x'VQ 



u = t-J2 + $^J3-$ r ^^dr, (77) 



where 



R = Wr = ^JQ, (78) 
e = G, N = H, T = U. (79) 

Jro VQ \rj Jro VQ Jro XVQ 



The case x = Transformations are clearly simplified when we choose the reg- 

ularizing factor X independent of the momenta G, H, and U. Thus, if x = xi'^)^ we 
have 

Xg = Xh = Xu = 0, (81) 

and the rightmost terms on the right side of (I7^ - (I77|) vanish, which turn these equations 
into 

/ = $fX3, (82) 
g = e + GIi + ^G^s, (83) 
h = iy + ^Hl3, (84) 
u = t-l2 + ^u^z- (85) 

Note that x = 1 x = ^"^ make I2 = I3 and X3 = — Xi, respectively. We plan to 
discuss these choices in detail below, focusing on a two-parametric family V = V(r; /i, b) 
of radial potentials. 
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4.2 Families of the isochronal potential 

In what follows we only deal with the transformation ( l82l) -( l85|) and, in addition, limit our 
study to the specific case of Henon's isochronal potential [211 E] 

V = (86) 

where b and /i are parameters. A relevant, particular case of the isochronal is the Keplerian 
potential V = —fi/r that occurs for 6 = 0. 

For the isochronal potential (1861) we find convenient to write (1731) as 

where a, p, and a, are certain functions of the momenta and parameters, which only will 
be determined after the regularizing parameter x has been chosen. Then, 

Q = ap(---)(--- 

\S Si/ \S2 S 



where Si > {s = \rr^ + W) > S2 are the two possible roots of the conic Q = 0: 

p a (1 — e^) , , o P 

si,2 = , = = a (1 ± e), = 1 - ^ < 1. 



4.2.1 Introducing auxiliary variables 

The roots Si^2, the extreme values of s, make natural the introduction of the auxiliary 
variables ip and 0, defined through 



y/r'^ + b'^ = a {1 — e cos ip) , dV r"^ + b'^ = ae sin ifj dip, (90) 

V;^^=— , df^=L=) =--sin<^d<^, (91) 
1 + e cos \ vr^ + V J P 



in order to ease the solution of quadratures in ( IHOl) . which we find convenient to write 

r + 6^ ^ / 1 



X. ^ r^df^). (92) 



X2 = / -j=d{V^^^), (93) 

Xa = -r'^d(^^\=r^d(Vr^). (94) 

For given values of the regularizing factor x, the quadratures above, (J92i) - (1M1) . may be 
integrated in closed form without need of specifying the new Hamiltonian, which remains 
as a formal function of the momenta $ = $(-F, G, H, U). 

Details on the families of canonical transformations generated by the cases x = 1, 
X = b + y/W+r^, and X = ^"^ ^ire given below. Besides, for the Keplerian case 6 = we 
recover Delaunay's [8], Levi-Civita's [251 EE], and Hill's [22] canonical transformations, 
as particular examples of our families. An elegant alternative for the derivation of these 
historic transformations may be found in [2]. 

14 



4.3 Case I: Delaunay's family of transformations x = !• 

From (1731) and ( IHTI) we obtain 



a = /X, 



Then, introducing the change of fl90l) in fl93l) . X2 is easily solved to give 



{ip — e sinip). 



(95) 



(96) 



Analogously, the change of (1911) is introduced in (1^ to solve the quadrature Xi. We 
obtain 

where we have introduced two new auxiliary variables 0i, 02, defined by means of the 
trigonometric relations p7] 



tan- 



tan — 
2 



1 + e — b/a /l — e 

A n— \ tan — 

^l-e-6/aVl + e 2 



1 + e — 6/a 
^ 1 - e - 6/a 2 ' 



1 + e + 6/a / 1 - e 1 + e + 6/a 

A n— \ tan — = A ——tan—. 

^l-e + 6/aVl + e 2 \l-e + b/a 2 



(98) 
(99) 



Substitution of Xi and X2 = X3, given by (!97|) and (p6!) . respectively, in (l82l) - (!85l) 
provides the family of canonical transformations of the isochronal potential to which we 
attach the name of Delaunay. 

Among the variety of possible choices of the new Hamiltonian $, a simplifying option 
is to take $ = [/ + \E'(F, G, if; /i, 6), that makes u to remain as the original time. In 
addition, the choice "ifn = 0, "^f = = A*"^ (—2^)^^'^, may be solved for to give 



$ = [/-- 



2 (F + G')2 

that maximally simplifies the transformation equations: 



/ 
9 

h 
u 



ip — esinip, 
6 



G 



2 ^/G^TWJ^ 2 



(100) 



(101) 

(102) 

(103) 
(104) 



while retaining the topology of the original Hamiltonian. Delaunay's selection = 0, 
= leads to 



<^ = U 



2F2 



(105) 
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that depends on fewer momenta and, therefore, constrains the topology of the original 
system to periodic solutions only — which may be adequate for a perturbation study like 

In the Keplerian case b = 0, V = —fJ'/r, the selection of the new Hamiltonian 
$ = F) does not constrain the range of solutions and the specific selection of (11051) 
provides the popular Delaunay transformation that, taking into account that 0i = 02 = 4>, 
is 

f = ip — esinip, g = 9 — (j), h = u, u = t, (106) 
where the most extended notation writes L = F, i = f . 

4.4 Case II: Levi-Civita's family of transformations x = ^ + \/&^ + 

Now, we write 

a = $ + /., a = ^, p=^ + 2h--. (107) 

2U <P + /i a 

The change (1911) is used to integrate Xi , (|^ , and the change (p7i) to integrate X2 and 
J3, dnSD and dMD. We get 

X. = , (108) 

2^ 2y^G2 + 46 (/x + $) 

X2 = ^(^-esin^), (109) 

X3 = , ^'^^ (110) 

y2f7 + 46 (/i + $) 



where 0i, 02, are the same auxiliary variables defined in (1951) and (|99i) respectively. Then, 
substitution of Xi, X2, X3, ( I108l) - (lll0p . in (I82l) - (l85l) provides the family of canonical trans- 
formations of the isochronal potential to which we tie the name of Levi-Civita. 
The, Keplerian case 5 = 0, 0i = 02 = 0, gives 

/ = "^^yi^^' (111) 

g = ^-0 + $G-L^, (112) 



h = z/ + $j,-Lv', (113) 



u + $ 1 
u = t- , ij) - esm^) + —^j). (114) 



The canonical transformation still remains undefined. We note in (11111) that f = ip when 
= V^2f7, a partial differential equation that may be solved to give 



$ = V2U 



F + C,iG,H,U) , (115) 
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where Ci is an arbitrary function of G,H, and U. The simple choice Ci = n provides the 
"first" Levi-Civita [25] transformatioiij 

f = ilj, g = e-(j), h = u, u = t- -^=={f -esinf). (116) 

Other possibihty is to force $(/ = (// + $)/(2f/) so that only periodic oscillations are 
introduced in the time scale (11141) . Then, 

<^ = -fi + VuC2{F,G,H). (117) 

with C2 arbitrary. The combination of both conditions (11151) and (I117P leads to 

<^ = V2u(f~-^ + C{G, H)) , (118) 



with C arbitrary. If we further take C = 0, we get 

F 

/ = ^, g = 0-(p, h = u, u = t + — e sin f. (119) 
that is the famous "second" Levi-Civita [26] transformation. 



4.5 Case III: Hill's family of transformations: x = 

We get now 

« = /i, a = ^^, p= + 26--. 120 

2U fj, a 

The required quadratures are solved with the same changes of variables as before, giving 

I = I - (]21) 
' ' 2VG^ - 2$ 2VG'^-2<l> + 46/i' ^ ' 

= ^^(^-esinV^), (122) 



which substitution in ( ]82l) -(l85l) provides the family of canonical transformations of the 
isochronal potential to which we assign the name of Hill. 

Focusing on the Keplerian case b = 0, (pi = (j)2 = (f), a. simplification choice that makes 
/ = is to select $ in such a way that ^p^i = —0, then 



2 



G^-{F + C)^\, (123) 

where C = C{G, H, U) is an arbitrary function. Then, 
$^ = -(F+C), <^g = G-{F+C)Cg, <^h = -{F'+C)Ch, <I>u = -{F+C) Cu. (124) 



^Note that Levi-Civita proceeds by scaling r by 1/V2U, or choosing x — r/V2U in our notation. 
This is equivalent to scaling our new Hamiltonian by the same factor, and yields to Levi-Civita's original 
Hamiltonian ^ = F while the transformation (|116p remain identical 
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The trivial choice C = makes g = 6, h = u, and u = t — ^[ip — e sin ip) / ^/SIP. 

Choosing C as a hnear combination of the momenta, simply ads to the angle variables. 
Specifically, when we set C = —G we obtairQ 



^ = F [G -]-F] (125) 



2 



and recover the TR-mapping [3ll [TT] 



/ = 0, g = e-f, h = iy, u = t--L={f-esmf). (126) 

V o(7'^ 



5 Conclusions 

A deep insight in the computation of canonical transformations is obtained through a 
general formulation of the Hamilton- Jacobi equation in the extended phase space that, 
besides, includes a regularizing function. As far as the transformed Hamiltonian may re- 
main formal, one can obtain families of canonical transformations instead of single sets of 
canonical variables. Then, particular transformations that meet the user's required char- 
acteristics are systematically computed from the solution of partial differential equations 
that may enjoy very simple solutions. 

Our formulation shows that different transformations in the literature belong to the 
same families. Furthermore, it reveals fundamental features of the transformations. Thus, 
among the worked examples, we identify families having as members known canonical 
changes of variables that originally were defined through implicit relations. 
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